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Abstract: We apply the formalism of geometrothermodynamics to the case of black 
holes with cosmological constant in four and higher dimensions. We use a thermody- 
namic metric which is invariant with respect to Legendre transformations and deter- 
mines the geometry of the space of equilibrium states. For all known black holes in 
higher dimensions, we show that the curvature scalar of the thermodynamic metric in 
all the cases is proportional to the heat capacity. As a consequence, phase transitions, 
which correspond to divergencies of the heat capacity, are represented geometrically as 
true curvature singularities. We interpret this as a further indication that the curvature 
of the thermodynamic metric is a measure of thermodynamic interaction. 
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1. Introduction 

One of the most interesting results of modern theoretical physics is its direct relation 
with many areas of mathematics like differential geometry. General relativity, for in- 
stance, can be considered mathematically as an application of differential geometry. 
Once a metric is given which is compatible with a torsion-free connection and satis- 
fies Einstein's equations, the corresponding curvature turns out to be a measure of 
gravitational interaction. This is a fascinating result that combines apparently differ- 
ent concepts of geometry and physics, allowing us to study gravity by measuring the 
curvature of spacetime. In fact, this result can conceptually be generalized to include 
all the field interactions that are known in nature. The electromagnetic, weak, and 
strong interactions can classically be described by using the Minkowski metric and a 
gauge connection. In all the cases, the resulting gauge curvature can be considered as 
a measure of the corresponding field interaction (see, for instance, [p]]). 

During the last few decades several attempts have been made in order to introduce 
differential geometric concepts in ordinary thermodynamics. Hermann || formulated 
the concept of thermodynamic phase space as a differential manifold with a natural 
contact structure. In the thermodynamic phase space there exists a special subspace 
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of thermodynamic equilibrium states. Weinhold [Q proposed an alternative approach 
in which in the space of equilibrium states a metric is introduced ad hoc as the Hes- 
sian of the internal energy. In an attempt to formulate the concept of thermodynamic 
length, Ruppeiner |3J introduced a metric which is conformally equivalent to Weinhold's 
metric. The study of the relation between the phase space and the metric structures 
of the space of equilibrium states led to the result that Weinhold's and Ruppeiner's 
thermodynamic metrics are not invariant under Legendre transformations H ||, i.e. 
the geometric properties of the space of equilibrium states are different when different 
thermodynamic potentials are used. This result clearly contradicts ordinary equilib- 
rium thermodynamics which is manifestly Legendre invariant. Moreover, the question 
whether the curvature of the space of equilibrium states can be considered as a mea- 
sure for thermodynamic interaction remained unanswered. This was particularly clear 
in the case of the thermodynamics of black holes where a flat thermodynamic metric 
can be transformed into a non-flat metric by means of a Legendre transformation |7]]. 

Recently, the formalism of geometrothermodynamics (GTD) was developed in order 
to unify in a consistent manner the geometric properties of the phase space and the 
space of equilibrium states |§. Legendre invariance plays an important role in this 
formalism. In particular, it allows us to derive Legendre invariant metrics for the space 
of equilibrium states. It has been shown that there exist thermodynamic metrics which 
correctly describe the thermodynamic behavior of the ideal gas and the van der Waals 
gas. In fact, for the ideal gas the curvature vanishes whereas for the van der Waals 
gas the curvature is non-zero and diverges only at those points where phase transitions 
take place. Moreover, in the case of black hole thermodynamics in four dimensions we 
have shown recently || that there exists a thermodynamic metric with non-vanishing 
curvature which correctly describes the thermodynamic properties of those black holes. 
The main goal of the present work is to show that for all known asymptotically anti-de 
Sitter black holes in all dimensions there exists a thermodynamic metric with non-zero 
curvature which correctly describes the structure of phase transitions as dictated by 
the corresponding heat capacity. Consequently, our main result is that the curvature 
of the space of equilibrium states can be used in a general manner as a measure of the 
thermodynamic interaction of black holes. 

The study of classical gravitational configurations on a background with cosmolog- 
ical constant has been intensified in the last few years. First, cosmo logical observations 
indicate that a positive cosmological constant could be responsible for the present ac- 
celeration of the Universe. On the other hand, a negative cosmological constant plays 
a distinguished role in the conjectured AdS/CFT correspondence, according to which 
5-dimensional solutions of Einstein equations with negative cosmological constant can 
be used to derive certain statements about quantum field theory in four dimensions. 
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In this context, the thermodynamic properties of black holes in an AdS background 
acquire especial importance since they give information about quantum field theory at 
non-zero temperature. Charged, rotating black holes in an AdS background are known 
explicitly only in four [[K]] and five dimensions [11]. Moreover, Reissner-Nordstrom-AdS 
and Kerr- AdS black holes are known in all dimensions [12|, O. The thermodynamics of 



these higher dimensional black holes has been a subject of intensive investigation due to 
its importance in the context of the AdS/CFT conjecture [14, 15, [l|J. In 2+1 gravity, 
the BTZ black hole presents an interesting phase transitions structure fl7 l. Charged 
topological AdS black holes and their phase transitions were analyzed in |fl8 |. One of 
the interesting results is that all the intrinsic parameters that characterize black holes 
in higher dimensions can be treated, with certain care, as thermodynamic variables of 
ordinary thermodynamics. 

This paper is organized as follows. In Section [| we review the general formalism of 
GTD for black holes in arbitrary dimensions and introduce a Legendre invariant metric 
in the thermodynamic phase space which is used to generate the geometric structure 
of the space of equilibrium states. In Sections |] and |5] we investigate the structure 
of the phase transitions of the Reissner-Nordstrom AdS and Kerr AdS black holes, 
respectively, and show the points where phase transitions occur are characterized by 
curvature singularities of the thermodynamic metric. In the final Section |7| we discuss 
our results. Throughout this paper we use units in which G = c = k B = h=l. 



2. Geometrothermodynamics 

In order to describe a thermodynamic system with n degrees of freedom, we consider in 
GTD the thermodynamic phase space which is defined mathematically as a Riemannian 
contact manifold (T, O, G), where T is a (2n + 1)— dimensional manifold, is a linear 
differential 1-form satisfying the condition A (dO) n ^ 0, and G is a non-degenerate, 
Legendre invariant metric on T. Here A represents the exterior product, d is the exterior 
derivative, and (dO) n = dQ A ... A dQ (n— times). The submanifold £ C T defined by 
means of a smooth embedding mapping ip : £ — > T such that the pullback ip*(Q) = 
is called the space of thermodynamic equilibrium states. A Riemannian structure g is 
induced naturally in £ by means of g = <p*(G). It is then expected in GTD H that 
the physical properties of a thermodynamic system in a state of equilibrium can be 
described in terms of the geometric properties of the corresponding space of equilibrium 
states £. 

To be more specific we introduce in the phase space T the coordinates Z A = 
(&,E a ,I a ) with A = 0, ...,2n, and a — 1, ...,n. In ordinary thermodynamics, $ corre- 
sponds to the thermodynamic potential, and E a and I a are the extensive and intensive 
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variables, respectively. The fundamental differential form can then be written in a 
canonical manner as G = c/$ — 5 ab I a dE b , where 5 ab is the Euclidean metric. The metric 
components in T can be in general arbitrary C 2 — functions of the coordinates, i. e., 
Gab = Gab{Z c ). This arbitrariness is restricted by the condition that G must be 
invariant with respect to Legendre transformations. This is a necessary condition for 
our description of thermodynamic systems to be independent of the thermodynamic 
potential. This implies that T must be a curved manifold || because the special case 
of a metric with vanishing curvature turns out to be non Legendre invariant. Although 
in general any n— dimensional subset of the set of coordinates Z A can be used to co- 
ordinatize the submanifold £, for the sake of simplicity we choose the subset E a as 
coordinates of £. Then the smooth mapping if : £ — > T is given in terms of co- 
ordinates as ip : {E a } i — ► Z A = {<5>(E a ),E a ,I a (E a )}. Consequently, the condition 
f*{Q) = can be written as the expressions 

d§ = 5 ab I a dE h 1 ^L = 6 ab I b , (2.1) 

which in ordinary thermodynamics correspond to the first law of thermodynamics and 
the conditions for thermodynamic equilibrium, respectively. We see that the specifi- 
cation of the mapping if includes the specification of the relationship $ = §(E a ) that 
is nothing more but the fundamental equation from which all the information about 
a thermodynamic system can be obtained |Tj|. The second law of thermodynamics is 
implemented in GTD as the convexity condition 

> . (2.2) 



8E a dE b 

To complete our construction we need a metric G. There is a large arbitrariness in 
the selection of this metric since at this level it is only demanded that it satisfies the 
condition of Legendre invariance. For the sake of simplicity we will use the following 
choice 

G = (d$ - 5 ab I a dE b ) 2 + (6 ab E a I b )( Vcd dE c dI d ) , Vab = diag(-l, 1, 1) , (2.3) 

where r] ab is a pseudo-Euclidean metric in £ . This metric is a slight modification of 
the metric G 11 presented in |J which was used there to generate the simplest Legendre 
invariant generalizations of Weinhold's and Ruppeiner's thermodynamic metrics. It is 
easy to show that the metric ( |2.3p is invariant with respect to the Legendre transfor- 
mation [^(J 

{$, E a , I a } — > {$, E a , I a } (2.4) 
$ = $ - 5 ab E a I b , E a = -I a , I a = E a , (2.5) 
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and when "projected" on £ by means of g = <p*(G) generates the thermodynamic 
metric 

^{ E w){^"e^ dE " iEd ) ' (2 6) 

Once the fundamental equation $ = §(E a ) is known for a given thermodynamic sys- 
tem, the explicit form of the thermodynamic metric g can easily be computed. If 
the curvature of the thermodynamic metric is to be considered as a measure of the 
thermodynamic interaction, the metric ( |2.(j| ) should be flat only for systems with no 
thermodynamic interaction. Moreover, phase transitions associated with divergencies 
of the thermodynamic interaction should correspond to curvature singularities. We 
will see that the metric Q2.6|) satisfies these conditions in the case of thermodynamic 
systems represented by black holes. 



3. Four dimensional Kerr-Newman-AdS black hole 

In the Einstein-Maxwell theory with cosmo logical constant A, which follows from the 
action 

S BM = 77j f dM-det^)] 1 / 2 (R - F^F^ - 2A) , (3.1) 
16vr J M i 

the most general solution representing a black hole configuration is given by the Kerr- 
Newman-AdS solution [plj that in Boyer-Lindquist-like coordinates can be expressed 



as 



A r / asin 2 6> , \ 2 A e sm 2 6 { , r 2 + a 2 \ 2 9 / dr 2 d6< 



ds d V +^^{ adt -— d *) + ' (a7 + a7 

(3.2) 

where 

/ r 2 \ a 2 
A r = (r 2 + a 2 ) M +-J - 2mr + q 2 , A e = 1 - — cos 2 9 , (3.3) 



p 2 = r 2 + a 2 cos 2 e , S = l-^. (3.4) 

This solution describes the gravitational field of a charged, rotating black hole with 
cosmo logical constant A = — 3 /I 2 , where I is the curvature radius of the AdS spacetime. 
The electromagnetic potential is given as 

a V r a aqr sin 2 6 

A = ~ ' 4? = — ^ — ( 3 - 5 ) 

with the angular, magnetic component A v that appears as a consequence of the rotation 
of the black hole. The horizons are determined by the (positive) roots of the equation 
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A r = 0. In particular, the outer horizon is situated at r = r + and corresponds to the 
largest root. 

The physical properties of this spacetime can be understood by considering the 
physical parameters entering the metric functions. The area of the horizon is a well- 
defined geometric parameter A = J y/g^geedtpdO that can easily be calculated at r = r + 
and yields A = 47r(r^_ + a 2 )/S. The surface gravity k can be derived, modulo a trivial 
additive constant, from the equation k a k b a — Kk b = evaluated at the horizon for a 
timelike Killing vector field k a . In the case under consideration we have that k a d a = 
dt+£ld v , where Q is the angular velocity measured by a non-rotating observer at infinity, 
so that the surface gravity is given as k = [3r\ + (a 2 + l 2 )r 2 + — (a 2 + q 2 )l 2 ]/[2l 2 r + (r 2 ^ + a 2 )}. 
The situation is more complicated in the case of the physical mass (total energy) M, 
angular momentum J, and electric charge Q, because these parameters are usually 
defined for asymptotically flat spacetimes. For asymptotically AdS spacetimes several 
definitions are possible and the issue has been clarified only recently by using the laws 



of black holes thermodynamics and the formalism of isolated horizons |T5|, [L6| . It turns 
out that it is necessary to measure the angular velocity with respect to an observer 
which is not rotating at infinity |14]]. Then the computation of the intrinsic physical 
parameters results in 

M=|, J=f, <J = i. (3.6) 

The connection to thermodynamics arises when one considers the Bekenstein- 
Hawking entropy in terms of the horizon area, i.e. S = A/4. It is then easy to 
derive the generalized Smarr formula for the KN-AdS black hole 

which is the fundamental thermodynamic equation. It relates the total energy M of the 
black hole with the extensive variables S, Q, and J. As in ordinary thermodynamics, 
in GTD it is the fundamental equation from which all the thermodynamic information 
can be derived. With the choice E a = {S,Q, J}, the corresponding intensive variables 
become I a = {T, <p, Q}, where is the electric potential and Q is the angular velocity. 
Furthermore, with this choice M corresponds to the thermodynamic potential. In 
this way, we have introduced all the coordinates Z A = {M, S,Q, J,T,<p,Q} of the 7- 
dimensional thermodynamic phase space T which, according to Eq. (|2.3|) , becomes a 
Riemannian manifold with metric 

G = (dM - TdS - 4>dQ - VldJ) 2 + (ST + (f>Q + QJ) (-dSdT + dQdcj) + dJdVl) . (3.8) 
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This is a non-degenerate metric with det(GAs) = (ST + QJ + (f)Q) 4 /16 and non-zero 
curvature. Moreover, the contact structure of T is generated by the fundamental form 
G = dM — TdS — cf)dQ — QdJ. At the level of the phase space T, the metric ( |3.8|) 
plays an auxiliary role in the sense that it generates a Legendre invariant metric for 
the space of equilibrium states S with coordinates {E a }. To this end, we introduce the 
smooth mapping 

ip : {S, Q, J} ^ {M(S, Q, J), S, Q, J, T(S, Q, J), <f>(S, Q, J), Q(S, Q, J)} (3.9) 

by using the fundamental equation ( |3T7j ) and the condition y?*(0) = so that on £ it 
corresponds to the first law of black hole thermodynamics dM = TdS + <pdQ + QdJ. 
This, in turn, can be used to compute the dual intensive variables corresponding to the 
temperature 



dM S 2 ( 1 7T fQ^\ ( 3_ 7T 7T 2 Q 2 \ IT J 2 

8S ~ 8Mvr 3 U 2 + S + S 2 \l 2 + S S 2 ) 2MS 2 



the electric potential, 



BM QSfl^ ^V (3 . u; 



dQ 2Mn \l 2 S S 2 



and angular velocity 

° " dJ ~ M \l 2 ' S 
Moreover, according to Eq.( |2.6| ), the metric structure of S is given as 



n 8M J (\ 7T . 

^ = ^T = T7 «+ a) ■ ( 3 - 12 ) 




M ss 











M QQ 


Mqj 





Mqj 


Mjj 



g = (SM S + QM Q + JMj) M QQ M QJ , (3.13) 



where subscripts represent partial derivative with respect to the corresponding coordi- 
nate. Notice that no cross terms of the form g$Q or gsj appear in this expression, which 
would be proportional to M S q or M S j, respectively. This is due to the special choice 
of the auxiliary metric ( |3.8| ). Indeed, the minus sign in front of the term dSdT in G 
leads to the disappearance of the cross terms of g = <p*(G) that involve the coordinate 
S, i. e., gsj and gsQ- Our choice of G is in agreement with the condition of Legendre 
invariance and was inspired by inspecting the expression of the scalar curvature R. 
In fact, R always contains the determinant of the metric g in the denominator and, 
therefore, the zeros of det(g) could lead to curvature singularities (if those zeros are not 
canceled by the zeros of the numerator). On the other hand, as we will show below, 
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the locations of the divergencies of the heat capacity coincide with the zeros of Mss- 
Then, the choice of the metric ( |3.8| ) has the purpose of generating a metric g whose 
determinant is proportional to Mss, leading to a one-to-one correspondence between 
the divergencies of the heat capacity and singularities of the scalar curvature. 

In the thermodynamics of black holes, phase transitions must play an important 
role. Due to the absence of a realistic, microscopic model for the entropy of black 
holes, a problem which is related to the absence of a theory of quantum gravity, there 
is no unanimity about the definition of phase transitions f2l| . Nevertheless, one can 



adopt the point of view of ordinary thermodynamics and search for singular points in 
the behavior of thermodynamic variables. Such an approach was realized by Davies 



22] , showing that divergencies in the heat capacity indicate the points where phase 
transitions occur. We follow Davies' approach in this work. From the fundamental 
equation ( |3.7| ) it is straightforward to compute the heat capacity for the KN-AdS black 
hole: 

c = = £L (3.i4) 

dT M ss 



S (j* + I + ZL ^~) (f + 



4tt 4 J 2 



J_ I vr , n 2 Q 2 \ ( 6_ i vr\ _ fn , 2n 2 Q 2 \ ( 3 , tt _ n 2 Q 2 \ , 8tt 3 ( -kJ^ _ qrp2\ 
S 2 J \l 2 ~+~ S) \S ~+~ S 2 J \l 2 S S 2 J ~+~ S 2 \ S 2 DJ ) 

Phase transitions are then determined by the roots of the denominator of C, i.e. Mss — 
0. 

On the other hand, for the curvature of the metric (13.131) to be a measure of the 



thermodynamic interaction in the KN-AdS black hole, it must reproduce the phase 
transitions structure dictated by the heat capacity Q3.14 ). To verify this property in an 



invariant manner we compute the scalar curvature of the thermodynamic metric ( |3.13|) , 
and notice that its denominator is given by 

D R = 4(SM S + QM Q + JMjf{M 2 QJ - M QQ MjjfM 2 ss , (3.15) 

whereas the numerator is a rather cumbersome expression that can not be written in a 
compact form. We see that the denominator is proportional to the determinant of the 
metric fl3.13p . At first sight, the singular points of the heat capacity that are situated at 
Mss — correspond to true curvature singularities where the volume element vanishes. 
However, this is valid only if the numerator of the scalar curvature does not eliminate 
the zeros of the denominator. A numerical analysis shows that in fact the singularities 
of the heat capacity coincide with the singularities of the scalar curvature. We first fix 
the value of the cosmological constant and the entropy in such a way that we get for 
the charge and angular momentum an interval where the heat capacity diverges. The 
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same process is then repeated for different combinations of values for the cosmological 
constant and entropy. As a result we find all the intervals where divergences occur. 
Around the divergent points of the heat capacity we then investigate the behavior of 
the thermodynamic scalar curvature and confirm that it becomes singular. We also 
noticed that the singular points coincide with the zeros of M$s so that in fact the 
curvature singularities are situated at the points where phase transitions take place. 
The characteristic behavior of the heat capacity and curvature is depicted in figures |l| 
and 0. 
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Figure 1: Characteristic behavior of the 
heat capacity C in terms of the angu- 
lar momentum J. The chosen values are 
A = -1, 5 = 1, and Q = 0.01. The di- 
vergencies indicate points of phase transi- 
tions. 



Figure 2: The thermodynamic scalar 
curvature R in terms of the angular mo- 
mentum J. The values of the remaining 
parameters are as in figure |]. The sin- 
gularity is located at the point of phase 
transition. 



In contrast to the above analysis, where the cosmological constant A has been 
treated as a fixed background parameter, it is possible to raise A at the level of an 
intrinsic parameter of the black hole and to consider it as an extensive thermodynamic 



variable . This is an interesting possibility that follows from the Kaluza-Klein reduc- 
tion of certain supergravity theories which are relevant in M-theory In GTD this 
possibility can easily be handled. In fact, in this case the phase space is 9-dimensional 
with coordinates Z A = {M, S, Q, J, A, T, 0, Q, L}, where L is the intensive coordinate 
dual to A. The construction of the Riemannian structure in T is straightforward, ac- 
cording to Eq. (|2.3| ). Furthermore, the coordinates of the space of equilibrium states 
can be chosen as E a = {S,Q, J, A} so that the corresponding thermodynamic metric 
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becomes 



g = (SM S + QM Q + JMj + AM A 



\ 



-M ss \ 
M QQ M QJ M qa 

Mqj Mjj M ja 

M QA M ja M aa J 



(3.16) 



The scalar curvature of this metric is again singular at the roots of Mss — so that, in 
principle, it can reproduce the structure of the phase transitions of the KN-AdS black 
hole. 



4. Reissner-Nordstrom-AdS black hole in arbitrary dimensions 



For a spacetime with arbitrary dimension D, the Einstein- Maxwell action with cosmo- 
logical constant can be written as 



S r 



1 



167T 



d D x[-det(g^)} 



M D 



1/2 



R - F^FV + 



(D-l)(D-2) 
I 2 



(4.1) 



where I is the characteristic length of the AdS background that determines the cosmo- 
logical constant by 

(D — 1)(D — 2) 



A 



2P 



(4.2) 



Then the metric for the RN-AdS black hole may be written in static coordinates as 



24, 25 



ds 2 = -f(r)dt 2 + f~\r)dr 2 + r 2 dQ 2 D _ 2 , (4.3) 

where dfl 2 D _ 2 is the metric on the unit (D — 2)— sphere. The function f(r) can be 
expressed as 

= + + (4.4) 



r D-3 1 r 2(D-3) 1 p ' 

where \x and q are intrinsic parameters related to the mass and charge of the black hole. 
This is an exact solution of Einstein-Maxwell equations with electromagnetic potential 



At — - 



D-2 



2(D-3)J r D 



1/2 



q 



(4.5) 



The outer horizon is situated at r = r + where r + is the largest root of the equation 
f(r) = 0. From this algebraic equation it follows that 



,D-1 



D-3 



+ 



(4.6) 
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Furthermore, the horizon area is given by 

.4 



D-2 



(4.7) 



where ui D _ 2 = 2ix (yD ~ 1 ^ 2 /T[(D — l)/2] is the volume of the unit (D — 2)— sphere. The 
calculation of the physical mass and charge can be carried out either by using an 
appropriate generalization of the Arnowitt-Deser-Misner (ADM), which includes the 
case of asymptotically anti-de Sitter spacetimes |26|, ^7]], or by using as a guide the laws 
of black hole thermodynamics [15] . The resulting quantities can be written as 

[D-2)u D _ 2 „ [2(D-2)(D-3)]^u L 



M 



Q 



Q 



16tt ^ 8tt 

After some algebraic manipulations which involve the expressions given above for hori- 
zon area in the form S = A/4, the mass, charge, and the parameter /i, we obtain 



M 



(g ~ 2K_ 
16tt 



AS 



D-l 
D-2 



1 



AS 



2 

D-2 



+ 



2vr 2 Q 2 



(D -2)(D -3)S 2 



(4.9) 



This is the fundamental equation for the RN-AdS black hole in arbitrary dimensions. 
In the special case D = 4 we recover the expression obtained in the last section. From 
the fundamental equation it is easy to derive all important thermodynamic variables. 
So we obtain the temperature 



T 



dM 
~dS 



47T \W 



AS 



(D 



+ (D-3) 



AS 



1^2 2vr 2 g 2 



(D-2)S\ 



(4.10) 



the electric potential, 



dM 



dQ (D - 3)S 



AS 



(4.11) 



and the heat capacity 



C ~ I dT 



(D-2)S 



D-l 



+ (D-3) 



is 



2 

D-2 



2n 2 Q 2 
(D-2)S 2 



M : 



D-l 



(D-3) 



is 



o- 2 + 2(2D-5)tt 2 Q 2 



(4.12) 



(D-2)S 2 



These thermodynamic variables are then considered as independent quantities at the 
level of the 5-dimensional thermodynamic phase space T which can be coordinatized 
by Z A = {M, S, Q, T, 0}. The Riemannian structure of T is determined in this case by 
the metric 



G = (dM - TdS - (j)dQf + (TS + 



-dTdS + d(pdQ) . 



(4.13) 
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The space of equilibrium states 8 can be chosen as being determined by the simple 
mapping ip : {S, Q} i— > {M(S, Q), S, Q, T(S, Q), 4>(S, Q)}, where the explicit depen- 
dence of the intensive variables is as given above. The thermodynamic metric on 8 can 
be computed by means of the pullback g = (p*(G) that yields 



g = (SM S + QM Q ) 



-M ss 

M QQ 



(4.14) 



As for the corresponding scalar curvature, we have performed an exhaustive analysis of 
the singularities and obtained a behavior similar to that of the heat capacity. Figures 
D and |] show a characteristic example of the singular behavior of the heat capacity 
and the thermodynamic curvature. Our analysis shows that the curvature singularities 
reproduce the structure of the phase transitions of the Reissner- Nordstrom black hole 
in arbitrary dimensions. 
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Figure 3: The heat capacity C as a func- 
tion of the charge Q for the 10-dimensional 
RN-AdS black hole. Here A = -1 and 
S = 100. Phase transitions are clearly 
seen as divergencies of C. 
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Figure 4: The scalar curvature R for the 
thermodynamic metric (4.14) as a func- 
tion of the charge Q. The values of the 
remaining parameters are as in figure ||. 
The locations of curvature singularities 
and phase transitions coincide. 



5. Kerr-AdS black hole in arbitrary dimensions 

In an arbitrary spacetime of dimension D, the Kerr-anti de Sitter metric is an exact 
solution to the equations R^ u + (D — l)l~ 2 g ^ u = 0, which in Boyer-Lindquist coordinates 
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can be expressed as 



13 



~2 i „2 



^ + 1 ^ *? -M i+ ^T{t "^'^ 2 ' (51) 



where N = [(D — l)/2] is the maximum number of independent rotational parameters 
aj in N independent, orthogonal 2-planes, m is a parameter related to the total mass 
of the black hole, the parameter e is defined as e = (D — 1) mod 2, t is the time 
coordinate, r is the radial coordinate, <pi are N different azimuthal angles, and /Zj are 
(N + e) direction cosines. Moreover, the functions appearing in the metric are 

N+e 2 N + e 2 N 

^E^.^E^ntr^ 2 ). (5.2) 

i=l i=l * j=l 

(2 \ JV 2 
i+ 72- n^ 2 +^)' s *= i -f- (5 - 3) 
^ i=i 

The outer horizon is situated at the radial distance r = r+, where r + is the largest 
positive root of the algebraic equation 

_e-2 



r. 



( v i +7Fjn( r + +a * 2 )- 2m=o ' ( 5 ' 4 ) 



which once solved can be used to find explicitly the value of the entropy S = A/A. As 
in the previous cases, the presence of the cosmological constant requires especial care 
for the determination of the physical parameters. Let us denote by M the physical 
mass of the black hole and by Jj the momentum corresponding to the angular velocity 
Qi, measured by a non- rotating observer at infinity. Using as a guide the laws of 



thermodynamics, it can be shown that [JL5 



4tt \ } \ ^— ' Ej 2 / ' Ar A 

■ %=\ • 

(5.5) 



,1=1 / \j=l J " ~r 1=1 



M-- 



N \ _1 N i N 2 



EiT> ^^J] 1 ^' forevenD, (5.6) 



vi=l / j=l J i=l 
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N ^ - 1 



Ji = ^°~ 2 \J[ Ejj , for arbitrary D . (5.7) 

The next step is the derivation of a Smarr-like formula which relates all the physical 
parameters entering the metric, i.e. M = M(S,Ji) or equivalently S = S(M,Ji). 
To this end it is necessary to express the outer radius r + in terms of the physical 
parameters. In case this is possible, we first select a representation for the construction 
of GTD. Let us consider the M— representation. Then, the extensive thermodynamic 
variables are S and J±. If we denote by T and fli the corresponding dual intensive 
thermodynamic variables, the coordinates of the phase space T can be chosen as Z A = 
{M, S, Ji,T,Qi} so that dim(T) = 3 + 2[(D — l)/2], when the maximum number of 
angular momenta is considered. The Riemannian structure of T is completed with the 
metric 

G= i^dM -TdS -^ftidJij + + ^ ^ J»J ^-dSdT + ^ d ^id Ji J . (5.8) 

The Riemannian submanifold of equilibrium states £ is defined by means of the smooth 
mapping ip : {S, Ji} i — > {M(S, Ji), S, Ji, T(S, Ji), £li(S, Ji)}, which induces the first law 
of thermodynamics dM = TdS + YliLi ^idJi and the (N + 1)— dimensional metric 



9= \ SM s + Y^JiMj}j 



f -Man \ 

(5.9) 

\ M JnJi ... M JnJn J 



Mj.j, ... M JiJn 




It is straightforward to see that the curvature scalar of this thermodynamic metric 
contains the term 



N ^ N 



SMs + ^JiMjA [detiM^fMls (5.10) 



i=l 



in its denominator, where My = d 2 M/(dJidJj), so that curvature singularities could 
appear at the points where the condition Mss = is satisfied. On the other hand, 
as mentioned above, in the mass representation the heat capacity can be expressed as 
C = Ms /Mss- It follows then that phase transitions that occur when the heat capacity 
diverges, could be represented as curvature singularities of the thermodynamic metric. 
This holds, of course, only if the expressions appearing in the numerator of the curvature 
do not eliminate the zeros of Mss- We note that in ( |5.10|) the term in round brackets 
cannot be zero because it can be shown to be proportional to the total mass M as a 
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result of Euler's identity. As for the determinant of My, its zeros, if any, can be found 
only if the fundamental equation M = M(S, J{) can be written explicitly. 

The determination of the fundamental equation turns out to be a non-trivial prob- 
lem. In fact, one way to determine it is to compute the solutions of the algebraic 
equation ( |5.4|) which, in general, is a polynomial of order D mr + . We were not able to 
find explicit solutions. However, the case of vanishing cosmological constant (I 2 — > oo) 
and only one non-vanishing rotational parameter, say, a± = a, can be manipulated 
explicitly. From Eqs.(pT3|)- (|5.7| ), we obtain for this specific case 

(r 2 + a 2 )r^ 5 - 2m = , (5.11) 

Af=-^(D/2-l), S = -f±(r 2 + + a 2 )r°-\ J — — ■ ( 5 - 12 ) 
A straightforward manipulation of these equations results in the expression 



D/2 — 1 f^ D _ 2 \ Q-3 / „ J 



7T 



r2 



M=^^(^j-^M+4vr^) . (5.13) 



D-2 



which constitutes the corresponding fundamental equation. In turn, the heat capacity 
can be expressed as 

(D - 2)S [3S 2 + 207T 2 J 2 - D(S 2 + 4tt 2 J 2 )] (g 2 + Air 2 J 2 ) 
~ 3S 4 + 2Ati 2 J 2 S 2 + 240vr 4 J 4 - D(5 4 + 48tt 4 J 4 ) ' ^ 5 ' 14) 

The corresponding thermodynamic metric of the space of equilibrium states reduces to 
the 2-dimensional metric 

g = (SM S + JMj) (-M ss dS 2 + MjjdJ 2 ) , (5.15) 

independently of the dimension D. The expression for the scalar curvature cannot be 
transformed into a compact form. Only the denominator can be shown to contain the 
expression 

(D - 3) [3S 4 + 24tt 2 J 2 S 2 + 240vr 4 J 4 - D(S 4 + 48tt 4 J 4 )] 2 , (5.16) 

which determines the phase transition structure of the heat capacity ( |5.14p . To see that 
the singular behavior of the scalar curvature coincides with that of the heat capacity we 
perform a detailed numerical analysis of both expressions. The characteristic singular 
behavior is depicted in figures [| and |[ For all analyzed regions a similar behavior was 
detected, showing that in fact the points where phase transitions occur are characterized 
by curvature singularities of the thermodynamic metric. 
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Figure 5: The heat capacity C as a 
function of the angular momentum J for 
the 10-dimensional Kerr black with one 
rotational parameter. Here A = and 
S = 100. 
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Figure 6: The scalar curvature R for the 



thermodynamic metric ( 5.15 ) as a func- 
tion of the angular momentum J. The re- 
maining parameters are chosen as in figure 



6. The role of statistical ensembles 



In the above applications of GTD of black holes, the starting point of the analysis is 
the fundamental thermodynamic equation. Nevertheless, it is known that the thermo- 
dynamic properties of black holes can drastically depend on the choice of statistical 
ensemble fl~4" , p4] , |25| , |32fl . The fact that different ensembles lead to different heat ca- 



pacities implies that the phase transitions structure depends on the statistical model 
under consideration. The question arises whether GTD is able to correctly handle this 
dependence. We will show in this section that the answer is in the affirmative. For 
the sake of concreteness and simplicity we will show this in the explicit example of the 
RN-AdS black hole in D = 4. 

Let us first consider the canonical ensemble which is characterized by a fixed charge. 
In this case, according to Eqs. (|4.9|) - (|4.11|) , the relevant thermodynamic variables can 
be expressed as 



2 V r\ P J ' 4tt \r + r 3 + P J ' 9 r + 



where we have used the relation S = 7cr^_. Then, the heat capacity Cq = (dM/dT)(_ 



- 16 - 



becomes 



Cr 



2nri 



3r 4 + + l 2 (r 2 + -Q 2 ) i 



(6.2) 



_3r\-l 2 (r\ -3Q 2 ) 

The analysis of this case in the context of GTD can be carried out as in Section 
[|. The metric G of the phase space T is given by Q4.13 ), whereas the metric g of the 
space of equilibrium states S coincides with Q4.14D and can be written explicitly as 



(SM S + QM Q ) {-M ss dS 2 + M QQ dQ 2 ) 



3r\ + l 2 (3Q 2 + r 



3r% + l 2 (3Q 2 



■ dS 2 + dJQT 



(6.3) 



This metric determines the geometric properties of the manifold £. In particular, the 
curvature scalar can be expressed in the compact form 



R 



48/M (r+, Q) 



[3r% + l 2 (r% + 3Q 2 )f[3r\ - l 2 (r 2 + - 3Q 2 ) 



(6.4) 



where 



Af{r+,Q) = 3rQ 4 (2/ 2 -21r 2 h )-3Z 2 r 2 h Q 2 (2r-5Z 2 r2+6r^)-r^(2r+3/ 2 r;-45r^). (6.5) 
The curvature singularities are then given by the zeros of the equation 

3r% -/ 2 (r 2 -3Q 2 ) = (6.6) 



which are also the blow-up points of the heat capacity ( |6.2j ). This shows that in the 
canonical ensemble approach the structure of the phase transitions coincides with the 
singular structure of the thermodynamic curvature. 

We now consider the grand canonical ensemble that corresponds to a fixed electric 
potential <fi. To obtain the grand canonical potential one can use the Euclidean action 



method as in |25[ or, equivalently, apply a Legendre transformation as in [|32| which 
interchanges the role of the variables Q and 0. That is to say, one introduces a new 
thermodynamic potential M by means of the Legendre transformation 



M = M 



(6.7) 



In [32|, the potential M is interpreted as the internal energy. Then, from Eq. (|6.1| ) we 
obtain 



On the other hand, calculating the total differential of M and using the first law 
of thermodynamics, we obtain dM = TdS — Qd<p so that the dual thermodynamic 
variables T = dM/dS and Q = —dM /dtp can be written as 



1 



1 



+ 



3r 2 



Q = <jrr+. 



47T7+ V ' ' P 

Furthermore, it is straightforward to calculate the heat capacity 

„ , ... , , , -3rl + lHl-^ 



(6.9) 







= 2nr 2 + 


[dT 



3r 2 -P(l- 



(6.10) 



Hence, the blow-up points that determine the phase transitions structure coincide with 
the solutions of the equation 



3rj 



l 2 (l 



. 



(6.11) 



We now analyze the case of the grand canonical ensemble in GTD. The new ther- 
modynamic potential is M = M(S, <p) so that the coordinates in the phase space T 
are Z A = {M, E a , I a } = {M, S, <j),T, —Q}. Then, the fundamental contact form can 
be written it its canonical form as = dM — TdS + Qdcf). The particular metric ( |2.3| ) 
we are using here for the geometry of the manifold T takes the form 



G = ( dM - TdS + Qdcj) ) -(ST- <f>Q) {dSdT + dcpdQ) . 



(6.12) 



The submanifold £ of equilibrium states with coordinates E a = {S, 0} is defined by 
means of the smooth map ip : £ — > T that in this case implies the explicit dependence 



if : {E a } i— > {Z A (E a )} = { M{S, 



),S,<f>,T 



dM n _dM 



(6.13) 



In turn, the geometric properties of £ are described by the metric g = <f*(G) which 
becomes 



g = [SM s + <j)Mt) [-M ss dS 2 + 



3r 2 + -l 2 (5(j) 2 - 1) 



3r 2 + -l 2 (l- 
8n 2 Pr 2 



2 i J2. Ji2 



dS z + r 



A straightforward calculation of the scalar curvature for this metric yields 

16/W(r+, 0) 



R 



r 2 + [3r 2 + - P(5(j) 2 - l)} 3 [3r 2 + - P(l - 2 )] 2 



(6.14) 



(6.15) 
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where 



A/"(r + ,0) = 5/ 4 4 (^ 2 + 24r2) + 3/ 2 r2 2 (6Or2 -47/ 2 ) + 3(/ 6 + / 4 r2 -6/ 2 r^ + 36r^) . (6.16) 

From the denominator of this expression and Eg. (|6 . 1 1|) , we conclude that there exist 
curvature singularities at those points where the heat capacity blows up. There is 
an additional term in the denominator which under the condition 3r+ — l 2 (5<f) 2 — 1) = 
can, in principle, lead to curvature singularities. However, it can easily be seen that 
this term is proportional to the conformal factor of the metric, i. e., SMs + 4>M^, which 
in turn, according to Euler's identity, is proportional to the internal energy M of the 
black hole. This seems to prohibit the appearance of additional singularities. 

We conclude that GTD can handle correctly the different thermodynamic schemes 
that follow from different statistical ensembles for the RN-AdS black hole in D = 4. 
Similar analysis can be performed for other types of black holes and we expect similar 
results. The fact that in general the denominator of the scalar curvature contains the 
denominator of the heat capacity is an indication that, even in the case of different 
statistical ensembles, there exists a correspondence between singular points of the heat 
capacity and curvature singularities. 

7. Discussion and conclusions 

The main result of this work is that in the space of equilibrium states of all known 
asymptotically anti-de Sitter black holes in arbitrary dimensions there exists a thermo- 
dynamic metric whose curvature is singular at those points where phase transitions of 
the heat capacity occur. This has been shown by considering a particular metric in the 
thermodynamic phase space, and applying the formalism of geometrothermodynamics. 
An important property of our choice of thermodynamic metric is that it is invariant 
with respect to Legendre transformations so that the properties of our geometric de- 
scription of thermodynamics are independent of the choice of thermodynamic potential 
and representation. 

The explicit examples considered in this work include the Reissner- Nordstrom black 
holes and the Kerr black holes on an anti-de Sitter background in arbitrary dimensions 
(D > 4). In the last example, we were able to explicitly analyze only the case of 
one rotational parameter with vanishing cosmological constant. The most general case 
of non- vanishing cosmological constant and arbitrary number of rotational parameters 
leads to a set of algebraic equations which we were not able to solve. Hence the 
fundamental thermodynamic equation cannot be written explicitly. Nevertheless, even 
in this general case our results show that if the numerator of the resulting expression for 
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the curvature scalar does not cancel the zeros of the denominator, then the curvature 
singularities are situated at those points where the heat capacity diverges, a fact that 
in ordinary thermodynamics is considered as an indication of a phase transition. The 
four-dimensional Kerr-Newman black hole on an anti-de Sitter background was also 
analyzed, obtaining similar results. The Kerr-Newman-AdS black hole is not known 
in higher dimensions. Thus, the examples considered in this work include all known 
higher dimensional solutions for black holes on an anti-de Sitter background. 

The starting point in our geometrothermodynamical analysis is the metric G of the 
thermodynamic phase space. In this work we used the particular metric ( [2.3|) which 
was chosen such that the metric g of the space of equilibrium states takes the specific 
form ( |2.6|) , whose determinant becomes proportional to <9 2 $ / '8E l dE l . In fact, the use 
of the metric rj a b, instead of 5 a b, in G guarantees that all non-diagonal terms of the form 
d 2 Q/dE 1 dE k , with k ^ 1, vanish. Furthermore, it is known that the scalar curvature 
always contains terms with the determinant of the metric in their denominator. It can 
therefore be expected that there exist true curvature singularities at those points where 
d 2 Q/dE 1 dE 1 = 0. We use in this work the M— representation, in which $ = M, and 
choose E 1 = S so that the singularities are expected at d 2 M/dS 2 = M ss = 0. On the 
other hand, in this representation the heat capacity can be written as C = Ms /Mss 
with divergencies at Mss — 0. For all known asymptotically anti-de Sitter black 
holes we have shown that the remaining terms of the thermodynamic scalar curvature 
do not cancel the zeros of Mss- Thus we conclude that the curvature of the space 
of thermodynamic equilibrium states can be interpreted in an invariant manner as 
a measure of the thermodynamic interaction. This is in contrast with other studies 
[ f2"g| , |2"m where the curvature of the space of equilibrium states depends on the choice 
of a specific representation so that, for instance, a flat thermodynamic metric can be 
associated to a system with non-trivial thermodynamic interaction. 

In the case of the RN-AdS black hole in D = 4, we analyzed the different thermo- 
dynamics which follow from the choice of different statistical ensembles. In particular, 
we considered the canonical ensemble, with fixed charge Q, and the grand canonical 
ensemble, with fixed electric potential <fi. We obtained the corresponding heat capaci- 
ties Cq and which lead to different phase transitions structures. According to our 
general procedure of GTD for black holes, we constructed explicitly for both ensembles 
the corresponding phase manifolds and the manifolds of equilibrium states. In both 
cases we obtained that the there exists a correspondence between singular points of 
the heat capacity and curvature singularities of the thermodynamic metric. We expect 
similar results in the analysis of more general black holes. The fact that in general 
the denominator of the scalar curvature contains the denominator of the heat capacity 
is an indication that GTD can handle correctly the different thermodynamic schemes 
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that follow from different statistical ensembles. 

Legendre invariance is an important element of our approach. It limits the number 
of metrics that can be used to describe ordinary thermodynamics in terms of geometric 
concepts. It is also essential in order to obtain results that are independent of the choice 
of extensive variables and thermodynamic potential. A different point of view, in which 
for a given thermodynamic system there exists a preferred thermodynamic potential 



30j pT|| , is necessary in order to explain the vanishing of Ruppeiner's thermodynamic 
curvature in cases where thermodynamic interaction is present as, for instance, in 
Reissner-Nordstrom black holes. We believe that ordinary thermodynamics, which is 
Legendre invariant, must not be changed when one tries to represent it in terms of 
geometric concepts. 

The computer algebra system REDUCE 3.8 was used for most of the calculations 
reported in this work. 
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